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Q Abstract 
O 

In the framework of QCD factorization, we study — > rfK + ^ decays. In order to more 

C ' 

reliably determine the phenomenological parameters Xh and Xa arising from end-point divergences 



in the hard spectator scattering and weak annihilation contributions, we use the global analysis 
for twelve B — * PP and VP decay modes, such as B — > irir, ttK, pir, pK, etc, but excluding the 
modes whose (dominant) internal quark-level process is 6 — > sss. Based on the global analysis, 



for twelve B — > PP and VP decay modes, such as B — » tttt, ttK, pir, pK, etc, but excluding the 

00 
OO 

■ we critically investigate possible magnitudes of Xh,a an d find that both large and small Xh,a 



(^T) , terms are allowed by the global fit. In the case of the large Xh a effects, the standard model (SM) 

o ■ 

prediction of the branching ratios (BRs) for B + ^ — > n'K + ^ is large and well consistent with the 
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D 



experimental results. In contrast, in the case of the small Xh a effects, the SM prediction for these 
BRs is smaller than the experimental data. Motivated by the recent Belle measurement of sin(20i) 



through B° — ► <f>K s , if we take into account possible new physics effects on the quark-level process 
b — ► sss, we can explicitly show that these large BRs can be understood even in the small Xh A 
case. Specifically, we present two new physics scenarios: R-parity violating SUSY and R-parity 
conserving SUSY. 
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I. INTRODUCTION 



From B factory experiments such as Belle and BaBar, copious experimental data on B 
decays start to provide new bounds on previously known observables with great precision 
as well as an opportunity to see very rare decay modes for the first time. There exist 
plenty of experimental data observed for charmless hadronic decays B — > PP (P denotes a 
pseudoscalar meson), such as B — > nn, irK, etc, and B — > VP (V denotes a vector meson), 
such as B — > pir, urn, pK, etc, which are well understood within the standard model (SM). 
However, among the B — ► PP decay modes, the BR of the decay modes B^ ^ — > rj'K^ ^ 
is found to be still larger than that expected within the SM. For last several years the 
experimental results of unexpectedly large branching ratios (BRs) for B — > rj'K decays have 
drawn a lot of theoretical attentions. The observed BRs for B ± — > rj'K in three different 
experiments are 

BBS 

B(B ± -> rj'K^) = (77.9l|^^) x 1(T 6 [BELLE], 

= (76.9 ± 3.5 ± 4.4) x KT 6 [BABAR], 

= (80^° ± 7) x 1(T 6 [CLEO]. (1) 

Many theoretical efforts have been made to explain the large BRs: for instance, ap pro aches 
using the anomalous g-g-rf coupling 0, 0, 0, Bl, high charm content in rj' lic|. the 

spectator hard scattering mechanism [UJ, ll2j, the QCD factorization (QCDF) approach 
13 1, the perturbative QCD (PQCD) approach [3| and approaches to invoke new physics 
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In earlier works on nonleptonic decays of B mesons, the factorization approximation, 
based on the color transparency argument, was usually assumed to estimate the hadronic 
matrix elements which are inevitably involved in theoretical calculations of the decay am- 
plitudes for these processes. This naive factorization approach ignores the nonfactorizable 
contributions from the soft interactions in the initial and final states. In order to compen- 
sate the nonfactorizable contributions, the naive factorization scheme has been generalized 
by introducing the effective number of colors iV c as a phenomenological parameter. In this 
generalized factorization, the renormalization scheme and scale dependence in the hadronic 
matrix elements has been resolved |20(. 

Theoretically, the QCDF approach has provided a novel method to study nonleptonic 
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B decays. In this approach, the naive factorization contributions become the leading term 
and as sub-leading contributions, radiative corrections from hard gluon exchange can be 
systematically calculated by using the perturbative QCD method in the heavy quark limit, 
where suppressed power corrections of 0(AQCD/ m b) are neglected. Since the nonfactorizable 
contributions in the naive factorization, such as the contributions from hard scattering with 
the spectator quark in the B meson and the contributions from weak annihilation, can 
be perturbatively computed, the phenomenological parameter N c used in the generalized 
factorization scheme is no longer needed to compensate the nonfactorizable contributions. 

However, in reality the b quark is not very heavy so that the power corrections in l/m&, 
particularly the chirally enhanced corrections, would not be negligible. The chirally en- 
hanced corrections come from twist-3 light cone distribution amplitudes (LCDAs), but un- 
fortunately the QCDF breaks down at twist-3 level because a logarithmic divergence appears 
in the hard spectator scattering at the end-point of the twist-3 LCDAs. A similar divergence 
also appears in the weak annihilation contributions. It is customary to phenomenologically 



treat these two end-point divergences by introducing model-dependent parameters |2l( : Xh 
for the hard spectator scattering contributions and Xa for the weak annihilation contribu- 
tions. Thus, it would be a less reliable case if these nonpertubative contributions of Xh 
and Xa become too large compared with the leading power radiative corrections. Since the 
prediction of the BRs for B —>■ PP and B — > VP decays strongly depend on the parameters 
Xh and Xa, it is essential to reliably estimate the effects of Xh and Xa- 

In this work we study the decay processes B^ ^ — > rfK^ ^ in the QCDF approach. In 
order to determine the parameters Xh and Xa more reliably, we use the global analysis 



as used in Ref. 



221 ] . However, our global analysis differs from that used in |22j], in the 



sense that we exclude the decay modes whose (dominant) internal quark-level process is 
b — > sss: for example, B — > <pK and B — > r/^M, where M denotes a light meson, such as 
7r, K, p, K*. The reason for excluding such modes is that the recent Belle measurement 
of the large negative value of sin(2</>i)0x s (0i is the angle of the unitarity triangle) through 
the time dependent decay process B° — > <pK s shows a possibility that there may be new 



physics effects on the quark-level process b — > sss |23J. Thus, to be conservative, in our 
global analysis within the SM, all the decay channels whose (dominant) quark-level process 
is b — > sss are excluded so that parameters Xh and Xa can be determined without new 
physics prejudice when using the global fit. For the analysis, we will use twelve B — > PP 
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and VP decay modes, including B — > nn, nK, pn, pK, utt, ujK. It turns out that both cases 
of the large and small X h ,a effects are allowed by the global fit. We will take into account 
both possibilities. In particular, motivated by the recent Belle result on sin(20i)0# s , we will 
seriously examine new physics effects on the large BRs for B — > rj'K. As specific examples of 
new physics models, we will present both R-parity violating (RPV) supersymmetry (SUSY) 
and R-parity conserving (RPC) SUSY scenario. 

This work is organized as follows. In Sec. II, we introduce the framework: the effective 
Hamiltonian for nonleptonic charmless B decays and the QCDF approach. The decay am- 
plitudes for B — > rj^K in the QCDF are presented in Sec. III. In Sec. IV, we discuss the 
global analysis for B — > PP and VP decays and calculate the BRs for B — > rj'K decays as 
well as B — > <f)K by using the inputs determined from the global analysis. We present the 
results for both cases of the large and small Xh,a effects. In particular, in the case of the 
small X h ,a effects, two new physics scenarios (RPV SUSY and RPC SUSY) are considered. 
We conclude the analysis in Sec. V. 

II. FRAMEWORK 

The effective Hamiltonian for hadronic charmless B decays can be written as 
(~i 10 

H*ff = -/={ E y P byp q [Ci(fj)o p M + c 2 (p)o p M + J2cmom 

V ^ p=u,c k=3 

-V tb V t * q [C 7l 7l + C 8g 8g ] } + H.c, (q = d, s) (2) 
where the dimension-6 local operators Oi are given by 



o\ = 
o u 2 = 
o 3 = 


(u a b a ) V - A (q~pUp) V - A , 

(u a bp)v-A(qi3U a )v-A, 

(qab a )v-A^2%q'p)v-A, 
q' 


01 = 

oi = 


(c a b a )v-A(qpCp)v-A, 

{c a bf})v-A{qpc a )v-A-, 

q' 


o 5 = 


{q a b a ) V -A ^{q'pq'^v+A, 

q' 


o 6 = 


(qpbah-A^i^'ph+A, 

q' 


7 = 


3 

q' 


8 = 


3 

-{qpb a ) V - A ^e q '{q! a q'p)v+A, 
q' 


o 9 = 


3 

-(q a b a ) V - A J2 e q'(QWp) v - A > 
q' 


Oi = 


3 

= ^(q(3b a ) V -Aj2 e q'^'a ( l'p)v~A, 

q' 


o 7l = 


g 


Os g = 


,^ mb q a a^(l +l5 )t a afS bpG% 



where q' denotes all the active quarks at the scale \i = 0(mb), i.e., q' = u,d,s,c,b. The 
operators Of, 2 are the tree operators, 0^-% are the strong penguin operators, Oj-io 
are the electroweak penguin operators, and 77 , 8g are the magnetic penguin operators. 
The Wilson coefficients (WCs) Cj(ju) are obtained by running the renormalization group 
equations from the weak scale down to scale p. We will use the WCs evaluated to the 
next-to-leading logarithmic order in the NDR scheme, as given in Ref. 24] . 

In the QCDF approach, in the heavy quark limit mj >> Aqcd ; the hadronic matrix 
element for B — » M\M 2 due to a particular operator Oj can be written in the form 



(MiM 2 |Oi|B) = (M 1 M 2 \O i \B) NF 



i+E'«(«.) B W— ) 



(4) 



where (M 1 M 2 |Oj|i?)NF denotes the naive factorization result. The second and third term 
in the square bracket represent the radiative corrections in a s and the power corrections in 
The decay amplitudes for B — > MiM 2 can be expressed as 



.4(5 -> MxM 2 ) = ^(S -> MxM 2 ) + A a {B -> AfiM 2 ) 



(5) 



where 



.u ; .u 2 ; 
AfiM 2 ) 



G p 
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V ^ p=u,c j=l 

^ fBfhhfM 2 E E^Vw & i 



(6) 
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E 

p=u,c i=l 

Here A*(B — ► MiM 2 ) includes vertex corrections, penguin corrections, and hard spectator 



scattering contributions which are absorbed into the QCD coefficients aj, and ^4 a (5 — > 
MiM 2 ) includes weak annihilation contributions which are absorbed into the parameter b^. 
For the explicit expressions of eij and 6j, we refer to Refs. j^l. I^. 

It is well known [2]] that both in the hard spectator scattering and in the annihilation 
contributions there appears logarithmic divergence in the end-point region. In Ref. 2jJ, 
Beneke et al. introduced phenomenological parameters for the end-point divergent integrals: 



X 



H,A 



r 1 dx / 

/ — =(l + p H: 

Jo X v 



1<PH,A 



hi 



m B 



(?) 



where Xjj and Xa denote the hard spectator scattering contribution and the annihilation 
contribution, respectively. Here the phases 4>h,a are arbitrary, 0° < 4>h,a < 360°, and the 



parameter pn,A < 1 an d the scale A h = 0.5 GeV assumed phenomenologically 



2l|. In 



principle, the parameters pn,A and <Ph,a for P — > PP decays can be different from those 
for B — > VP decays. Thus, for B — > PP and VP decays, from the end-point divergent 
integrals, eight new parameters are introduced: Ph^ai 4>W,a f° r B ~^ PPi an d Ph P a-i ^h^a fo r 
B -> VP. 



III. DECAY PROCESSES B±(°) -► r/^K ±( -^ IN THE QCDF APPROACH 



The decay amplitudes for P~ -> ^W-K" - and P° -> ?? (/) P; in the QCDF are given by 



■[V ub V: s (a[ + or + <o' + K ' + < W 
+K 6 K;(«4 ' + a?o + K ' + «s ')^i)l 



ru 
J V (>) 



\ ;,!,] ;„ I a 2 + 2a 3 - 2a 5 - ^(a 7 - a 9 ) - - i<£ ) P; 



+V/ c5 K; ^2a 3 - 2a 5 ~\{a 7 - a 9 ) - - ^ P 3 
+/*<o KfeV;* s (a 3 + < - a 5 + -(a 7 - a 9 - ai ) + - T^s) ^3) 
+ KbKt ( a 3 + a c 4 - a 5 + ^(a 7 - a 9 - ai ) + ^ajj - ^ajj^ P 3 



G 



— 1 



^JbJk (#„ + #) [WA + ( W£ + v cb v; s ) (h + b 

V2 



3 J 



,(8) 



A{B»^rPK°)= - i^fKF B ^ (,) (m 2 K )(m 2 B -ml (/) ) 



\v ub v:M' +K' -\< W 



{/Jo [K 6 C (02 + 2a 3 - 2a 5 - ^(a 7 - a 9 ) - - ±a$) P; 



CLr — —ds 1 P 3 



+Vd,V* (2a 3 - 2a 5 - ^(a 7 - a 9 ) - (a* - - 
+/*(0 V ub V* s (a 3 + al - a 5 + ^(a 7 - a 9 - ai ) + (< - T^s) ^3) 
+ V cb V* s (a 3 + a%- a 5 + ^(a 7 - a 9 - a 10 ) + (a£ - ^ajj^ P 3 ^) J 



6 



■^§/b/* too + /Jo) (Kk^ + WE) ( & 3 - i« 



} 3 



where 



R 



1(2) 



The coefficients of and 6j are expressed as 

C 2 



2m 2 K(0) 



(m b - m u{d) )(m u{d) + m s 



R* 



2m 2 , 



2m s (m b - m s ) 



a? 



«2 



C 2 + 



N c 

as = C i + - 



l + ^(v M + ^H(BM u M 2 ) 

47T V JS r 



a 5 



«5 W 



_p (') 
a 8 



C 4 + 
C 5 + ^ 

c 6 + — 

a 7 = C 7 + 
C 8 + 



C3 

N c 
Ch 
N c 
Cs 
N c 
Cr 



1 + 

1 + 

1 + 

1 - 



C F a s 

An 
C F a s 



V M + ^H(BM U M 2 ) 



. .V M + ^H(BM 1 ,M 2 ) 



1-6 
1 - 



C F Ot s 

An 

C F a s 

An 

C F a s 



V M + ^-H(BM x ,M 2 ) 



+ 



C F a s pP 

Attn/ m - 2 ' 



Vm + 12 



4tt z 



H(BM 1 ,M 2 



Cpo^ 
+ AnN c M ^ 
An' 



a 9 



a io 



r 4. 6,10 



c 



10 



Cg 
iV, 



1-6 
1 + 
1 + 



An 

%^ I + 12 + ^—H(BM U M 2) 
An \ N c 

C F a s 



An 

C F a s 

An 
C F a s 
An 



+ 



a 



9nN r 



-}P,ew 
M,3 ' 



V 



An 



M 



■H(BM 1 ,M 2 



+ 



a 



9nN r 



-P 



p,ew 
M,2 ) 



6 2 
63 



°3 



2^1 , 



cv 

^[C 3 ^ + ^(C 5 + iV c C 6 )], 

c 

(7F [C 9 ^ + ^(C 7 + iV c C 8 )], 



(9) 



(10) 



(11) 



N 2 

c 

where the superscript p is u or c, and the color factor C F = (N 2 - 1)/(2N C ) with N c = 3. 
The vertex parameter Vm and the hard spectator scattering parameter H(BM\, M 2 ), and 
the weak annihilation parameters A 1 , A* are given by 



V M = 12 In— -18+ / dx g(x)$ M (x), 
H Jo 



7 



H(BM 1 ,M 2 ) 



»i 

x / dy 







o € Jo (1 — x) 



(i-y) 



.4' 



7T 
1 



m 6 

2\ 



X (1 
2 



y) 



18 ( X A - 4 + — ) + 2r z x X A 



12na s r x X A {2X A - 1), 



(12) 



where g(x) 



l-2x 



\nx — 3«7r and fip 



{nil an d m 2 are current quark masses 



X—x J """ f~j- mi _|_ m2 

of the valence quarks of the meson P) and the chirally enhanced factor r y 



For the 



,K 



r x as in Ref. 



24 1 . 



1 — jf 

X A = Jg 1 — is a logarithmically divergent integral. For the wave function $b(£) of the B 
meson, we take the following parametrization: 



o 



(13) 



where the parameter A# is estimated as Xb = (350 ± 150) MeV 2l|. For the K and rf 
meson, we use the asymptotic forms of the LCDAs (2! ]: 



$,,/ = 6x(l — x), 



(14) 



where $m(^) and $^ f (x) are the leading twist LCDAs and twist-3 LCDAs of the meson 
M = K, T)', respectively. The explicit expressions of the QCD pen guin parameters Pm% an d 



the electroweak penguin parameters Pff™ can be found in Refs. |2l|, |24j. The coefficients 
Oj and a[ in Eqs. (jHJ) and Q include the different vertex and hard spectator scattering 
contributions: for a i; Vm = V v > and H(BM\, M 2 ) = H(BK,i]'), while for a^, Vm = V# and 
H(BM h M 2 ) = H(Br)',K). 

Note that in Eq. (fT2*|) the hard spectator scattering parameter H(BM\, M 2 ) includes a 
logarithmically divergent integral Jq dy/(l — y) which arises from the twist-3 contribution, 
and the weak annihilation parameters A 1 and A? include another logarithmically divergent 
integral X A . 

For the rj — 7/ mixing, we use the following relation: 



\r)) = cos 9 s \r) S ) - sin o \r) O ), 
|t/) = sin 9 8 \rj 8 ) + cos 9 \rj ), 



(15) 



S 



where r/s and r] are the flavor SU(3) octet and single, respectively. The mixing angles are 
6*8 « —22.2° and 9q ~ —9.1° 25]. The decay constants and form factors relevant for the 



B — > i]^ transitions are given by 



_ cos98 __ sm9o , /< = _ 2 „„,«, 



= fff -5^) , <?' = f~ + ) . (16) 

IV. GLOBAL ANALYSIS AND NUMERICAL RESULT 

In order to calculate the BRs for B decays in the QCDF approach, various input pa- 
rameters are needed, such as the CKM matrix elements, decay constants, transition form 
factors, LCDAs, and so on. Among those input parameters, it is urgently essential to reliably 
estimate the annihilation parameter Xa and the hard spectator scattering parameter Xh'- 
more specifically, pa, (ft a, Ph, and (fin, because the predicted BRs strongly depend on the 
parameters Xa and Xh (see Figs. 1 and 2). Unfortunately, within the QCDF scheme, Xa 
and Xh are purely phenomenological parameters, so there is no definite way to determine 
them. Therefore, in order to determine the values of pa^^a, Ph, and more reliably, 



in this work we follow the global analysis, used in Ref. 



22 1 . For the detailed discussion 



on the method of the global fit, we refer to Ref. |22j]. As explained in Sec. I, different 
from the global analysis used in j^], we do not include the decay modes, such as B — > <pK 
and B -> r/ (/) M (M denotes a light meson: e.g., tt, K, p, K*), whose (dominant) internal 
quark- level process is b — > sss. In this way, within the SM, the parameters pa, (ft a, Ph, (fin 
can be determined without new physics prejudice when using the global fit. Specifically, we 
use twelve decay modes, such as B — > nn, irK, pir, pK, un, and uK, as listed in Table I 

□ □□ 

First, we examine the dependence of the BR for B + — > rj'K + on the effects of Xa and Xh- 
Figure 1 shows B(B + — ► r)'K + ) versus $ P A P (solid line) or <p H p (dotted line). In each case, 
4>a P or $h P varies from to 2tt. For the solid line, other inputs are set as p p A F = p H P — 1> 
<P H P = -23°. For the dotted line, p p A p = p p H p = 1, <p pp = 57°. We see that the predicted BR 
for B + — > r]'K + strongly depends on (p PP and <Ph P - In particular, as the value of § P A P varies, 
the predicted BR can change by a factor of about 2.5 (e.g., from 45 x 10~ 6 to 116 x 10~ 6 ). 
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FIG. 1: Dependence of the BR for B + — » 7]'K + on ^ pp (solid line) or 0^ (dotted line). Here the 
following values of the other parameters are used: p pp = p H p = 1 (for both lines), 4> H P = —23° (for 
the solid line), PP = 57° (for the dotted line). The shaded region is allowed by the experimental 
data. 

The allowed values of PP are in certain narrow regions which can be practically found by 
the global analysis. Similarly, Figure 2 shows B(B + — ► i]'K + ) versus p pp (solid line) or 
p^f (dotted line). In each case, p pp or p^f varies from to 1. The other inputs are put 
as PP = —23° and <f> H p = 57° for both lines, p^ = 1 for the solid line, and p pp = 1 
for the dotted line. The predicted BR for B + — > r]'K + is also dependent on p pp and p^f ', 
but its dependence on p PJ ^ is weaker than that on PP . We notice that the prediction of 
B{B + — > r]'K + ) is very sensitive to the effect of Xa through PP . This feature also holds 
for the neutral mode B° — ► r]'K°. 

We find that the best fit [and also the "good" fit (see the discussions below [Case 1])] 
of the global analysis favors large effects of the parameters Xa and Xh- This tendency 



22, m. 



is consistent with the results of other previous works done in the QCDF scheme 
But, as mentioned in Sec. I, if the nonperturbative effects of Xa and Xh are too large or 
dominant compared with the leading power radiative corrections, the theoretical predictions 
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FIG. 2: Dependence of the BR for B + — ► rj'K + on (solid line) or (dotted line). Here the 
following values of the other parameters are used: 4>^ p = —23° and 0^ = 57° (for both lines), 

experimental data. 



' f — 1 (for the solid line), p^ P = 1 (for the dotted line). The shaded region is allowed by the 



based on these effects would be less reliable and become questionable. Therefore, one can 
seriously ask the following question: Is it possible to find a global fit where the effects of 
Xa and Xh are rather small (so the theoretical predictions based on these effects would be 
more reliable), but its \ 2 value is still acceptably small? In fact, it turns out that such an 
acceptable fit with the small effects of Xa and Xh can be found. 

For calculation of the BRs for — > r]'K +( -°\ we take into account two different 
possibilities as discussed above: [Case 1] with the large Xa and Xh effects (favored by the 
best and "good" fit, but less reliable), [Case 2] with the small Xa and Xh effects (more 
reliable). 

[Case 1] with the large Xa and X H effects 

We first try to find the best fit of the global analysis for the twelve B — ► PP and VP 
decay channels shown in Table I. Our result shows that based on the theoretical inputs for 
the best fit (with xLin = 7.5), the predicted BR for B + — > r/'K + is consistent with the 
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TABLE I: Experimental data and the "good" fit values of CP-averaged branching ratios (in unit 
of 10~ 6 ) for B — > PP and VP decays used in our global analysis. 



Decay mode 


Weighted average(Exp.j 


lit 


Decay mode 


Weighted average(Exp.J 


1 it 


D — > 7T 7T 


K /I9 -|- n 8^ 




o0 , _+_- 

-D — > 7T 7T 


A Z.K _l_ n A A 


A 7^ 




20.8 ± 1.4 


20.5 




18.1 ±0.8 


18.9 


-► TT°K+ 


12.7 ± 1.1 


11.4 


B° vr°K 


11.2 ± 1.4 


8.5 


B+ - 7r+p° 


8.6 ±2.0 


7.6 


5° - tt V 


22.7 ±2.5 


23.5 


S+ -► Ll>7T+ 


6.0 ±0.9 


6.45 


5° -► K+p- 


8.0 ±1.7 


9.5 


5 + -► wK+ 


5.6 ±0.9 


5.25 




5.3 ± 1.5 


4.45 



experimental data, but the prediction of B(B + — > <f)K + ) is too small compared with the 
data (for the data, see Table II), as the best fit predicted BRs are 

B(B + -> r]'K + ) = 74.7 x lCT 6 , 

B(B + -> 0fT+) = 4.02 x 10~ 6 . (17) 

It happens because the internal interference between different contributions (e.g., contri- 
butions from the hard spectator scattering and weak annihilation) to the decay amplitude 
for B + — > i]'K + is quite different from that for B + — > <f)K + (for example, see Table II). It 
turns out that it is possible to obtain successful fits to all B — > rj'K and B — > <f>K data, if 
one assumes that there are new physics effects on the quark-level process b — > sss. We will 
discuss this possibility later in SUSY scenarios. 

Since the input parameter values for the best fit are not consistent with the experimental 
result such as the BR for B — > <f>K, we investigate another possibility that there may exist a 
"good" fit for which the predictions based on the inputs are consistent with the experimental 
measurements including B — * rj'K and B — > <f)K, and whose Xmin value is still quite small. 
In fact, we find such a "good" fit with xLin = 8-6 for the twelve decay modes. Notice that 
this Xmin value is not much different from that of the best fit. In Table I, we list the "good" 
fit values of the BRs for the relevant B — > PP and VP decay modes. The corresponding 
theoretical inputs are given by 

A = 0.2205, A = 0.814, 3 = 72°, \V ub \ = 3.49 x 10" 3 , 

H = 2.1 GeV, m s (2 GeV) = 85 MeV, f B = 220 MeV, X B = 200 MeV, 
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TABLE II: Experimental data and the prediction of the branching ratios (in unit of 10 -6 ) for 
B — ► -q'K and B — ► 0.K" decays. Here the inputs for the "good" fit are used. For comparison, the 
predicted BRs for three cases are also listed: (i) for Xa = Xh = 0, (ii) for only Xa = 0, (iii) for 
only Xh = 0. 



Decay 


mode 


Exp. data 


Prediction 


X A 


= X H = 


Xa = only 


Xh = only 


B+ -» 




77.6 ±4.6 


78.5 




52.3 


64.9 


66.6 


B ^ 




65.0 ±6.0 


71.6 




47.8 


59.5 


60.7 


B+ -► 


cj)K + 


9.3 ±0.7 


8.85 




2.27 


1.51 


6.49 


5°^ 


cj)K 


8.2 ± 1.1 


8.01 




2.06 


1.37 


5.85 



F Bw = 0.23, R nK = 1, A Bp = 0.31, 



pY = pY = pY = pV = i, 

07 = 57°, 07 = 52°, 0^ = -23°, 0^ = 180°, (18) 

where A = |Ks| and R nK = {UF BK )/{f K F m ). The parameter A is defined by AX 2 = \V cb \ 
and 03 is the angle of the unitarity triangle. J'b is the B meson decay constant, and F 3 ™ 
and A Bp are the form factors for the transition B — > 7r and i? — > p, respectively. 

Indeed the BRs for -> r/iT^ ) and -> 0iT f (°) calculated by using the above 

inputs are in good agreement with the experimental measurements as shown in Table II. 
Therefore, our result shows that the large BRs for the processes B + ^ — > i]'K +( - ^ as well as 
the BRs for -»■ 0AT+(°) can be consistently understood, based on the global analysis 

for B — > PP and decays, where the values of the pure phenomenological parameters 
pYh-i pYh-i ^Yh and <PYh are reasonably determined. The BRs for B + ^ — > ^i^" 1 "^ ^ are 
estimated as (1 ~ 2) x 10~ 6 which are also consistent with the data ODD B(B + -> 
?7i^ + ) = (3.7 ± 0.7) x lO" 6 and B(B° -> t^FT ) = (2.9 ± 1.0) x 10" 6 . 

However, we note that the inputs given in Eq. ()18|) provide large effects of Xa and Xjj: 
e.g., p BP H = pYh = 1 [ see Eq. ©]• In order to explicitly estimate the effects of Xa and 
Xh, we also examine three interesting cases. In the fourth column of Table II, the BRs for 
B — > vjK and B — > <pK are calculated for Xa = Xh = 0. Similarly, in the fifth and last 
column, those BRs are calculated under the assumption of Xa = or Xh = 0, respectively. 
We see that the contributions from the terms involving Xa and Xh are quite large for both 
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B — > rj'K and B — > 0if decays. In particular, for B(B +< ^> — ► the contribution of 

(i.e., weak annihilation contribution) dominates over all the other contributions. /,From 
the table it is clear that the internal interference between the effects of Xa and Xh on 
B — > rj'K is constructive, while that on B — ► is destructive. It should be stressed that 
in this scenario (i.e., with the large effect of Xa,h allowed by the "good" fit), there is no 
room for invoking new physics effects on the quark- level process b — > sss, which is implied 
by the large negative value of sin(20i) recently measured by Belle 2^. 
[Case 2] with the small Xa and Xh effects 

As already emphasized, if the nonpertubative contributions of Xa and X H are too large, 
the predictions based on these contributions become less reliable and suspicious. However, in 
[Case 1], we noticed that the contribution of Xa is very large, especially for B + ^ — > <j)K + ^ 
modes. Therefore, it is natural to investigate presumably more reliable scenarios, where the 
effects of Xa and Xh are rather small or at least not dominant. 

Using the global analysis for the twelve decay modes shown in Table I, we find such a 
fit (with xLin = 18.3) with the (relatively) small Xa and X H effects. The corresponding 
theoretical inputs for this fit are as follows: 

A = 0.2198, A = 0.868, <p 3 = 86.8°, \V ub \ = 3.35 x 10~ 3 , 
// = 2.lGeV, m s (2 GeV) = 85 MeV, fg = 220 MeV, 
F Bn = 0.249, R nK = 1, A Bp = 0.31, 
P 7 = 0, PY = 0.5, /&> = !, pY = 0-746, 

<PY = -6°, 0^ = 0^ = 180°. (19) 

Note that in this case the effect of the weak annihilation parameter Xa is relatively small 
(i.e., Pa P = and p v A p = 0.5), and the effect of the hard spectator scattering parameter Xh 
is very small, because p H P = 1, pif = 0.746, and <ft H P — <Ph P = 180° so that the terms 1 
and ph^ h in Xh [see Eq. (J7J)] cancel each other. 

Based on the above inputs, the BRs for B — > rj'K and B — > <pK are predicted as 

B(B + -> t]'K + ) = 51.1 x 10" 6 , B(B° -> rj'K ) = 46.8 x 10" 6 , 

B(B + -> 0X+) = 7.29 x 10" 6 , B(B° -> 0ir o ) = 6.65 x 10" 6 . (20) 

These BRs are quite small, especially for B — > rj'K, compared with the experimental data, 
because of the small effects of Xa and X H as well as the other fitted parameters such as </> 3 . 
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Since both processes B — > rj'K and B — > cf)K have the same (dominant) internal quark-level 
process h — ► sss, we take into account the possibility that there could be new physics effects 



on the process b — > sss: for instance, as considered in [17|, 2t| in order to explain the large 
negative value of sm(2(pi) ( j ) K a reported by Belle. We investigate whether it is possible to 
understand the difference between the BRs given in Eq. (J2U)) and the experimental data, by 
invoking new physics. 

As specific examples, we consider two new physics scenarios: R-parity violating (RPV) 
SUSY and R-parity conserving (RPC) SUSY, 
(a) R-parity violating SUSY case 

The RPV part of the superpotential of the minimal supersymmetric standard model can 
contain terms of the form 



(21) 



where Ei, Ui and Di are respectively the i-th type of lepton, up-quark and down-quark 
singlet superfields, Lj and Qi are the SU(2) L doublet lepton and quark superfields, and H 2 
is the Higgs doublet with the appropriate hypercharge. 

For our purpose, we will assume only /'—type couplings to be present. Then, the effective 
Hamiltonian for charmless hadronic B decay can be written as jl5j], 



H*f f (b -> djd k d n ) 

Heff(b -> UjU k d n ) 



d R 

jkn 



U jkn 



dnalL^jp dkp'-y^nba 

UkcflL U jP dnPl^Rbc 



4- H L 

' u jkn 



(22) 



Here the coefficients d L ^ k R and uf kn are defined as 



d R 

u jkn 



U jkn 



3 // //* 

\ - L ijk L in?> 

3 ;/ //* 

^ l ijn L ik?> 



d L 

u jkn 



3 // //* 
\ ~> i i3k i inj 

k s<7 



{j,k,n = 1,2) 



(j,k = l, n = 2) 



(23) 



where a and (3 are color indices and 7^ L = 7 M (1 ± 75). The leading order QCD correction 



to t 
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lis operator is given by a scaling factor f — 2 for = 200 GeV. We refer to Refs. 



16( for the relevant notations. 



The RPV SUSY part (relevant to the quark-level process b — > sss) of the decay amplitude 
of B~ — > rj'K" is given by 



/.RPV 



(d 



L -d R 

222 u 222 



a 4 — a 5 ) + A^,a,i 



(24) 
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where 



mz 



m = 



Here the coefficients are expressed as 



m. 



A { 
v' 



u(s) 



j<s) F B^K { 



(25) 



a i 



CfCKs 

AnN r 



1 1^-0,(0) 

3 jj, 



°5 = TT 



1 

N, 



C F a s 



4vr 2 

12 + — H(BK, rf) 



«6 



«6 



1 + 



^-6,(0) 

3 /i 



-In Gjc(O) 

3 



(26) 



16 I 27T,' 

9 ^ 3 l - 



where G x (0) = | + f z and G x (0) 

It has been noticed |l7J that A^F"^ is proportional to ((i 222 — ^222)' while the RPV part 
of the decay amplitude of B — > <pK is proportional to (c? 222 + ^222)- ^ nas been also pointed 
out 17] that the opposite relative sign between g? 222 an d ^222 m the modes B — > rj'K and 
B — > <pK appears due to the different parity in the final state mesons rf and <p, and this 
different combination of (c? 222 — c^ 22 ) and (<i 222 + ^222) m these modes plays an important 
role to explain both the large BRs for B — > rj'K and the large negative value of sin^^),^ 
at the same time. 

We define the new coupling terms <i 222 and df 22 as follows: 



^222 c*- lA'^aA'* 



v i32 /v i22l 



d 



222 



K l^22 A i23l e 



i6 h 



(27) 



where 6l and Or denote new weak phases of the product of new couplings A^ 32 A^ 22 and 



•M22-M:23> respectively, as defined by A 3 3 2 A 3 * 22 



1^332^322 I e dIiu A 322 A 323 



A009 A 



I* 



322^3231 



r]'K + ^ and B+<® 



We find that the experimental measurements of the BRs for B + ^ 
(f>K + ^ can be consistently understood for the following values of the parameters: 



|A 322 | = 0.076 , I A^l = 0.076 , |A 323 | 



0.064 



6 L = 1.32 , B R = -1.29 , m S usY = 200 GeV. 



(28) 



Our results are summarized in Table III. In addition to the parameters given in Eq. (|2*K|) . we 
also used the additional strong phase 5' = 30°, which can arise from the power contributions 
of AgcD/nT'b neglected in the QCDF scheme, and whose size can be in principle comparable 
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TABLE III: The branching ratios (in unit of 10~ 6 ) for B — > rj'K and B — ► 0if decays calculated 
in the framework of R-parity violating SUSY. Here the inputs for the fit with small Xa and Xh 
are used. 



Decay 


mode 


Prediction 


Decay 


mode 


Prediction 


B+ 




74.0 


B° — > 


t/K° 


67.7 


-> 




10.2 


B ^ 




9.5 



to the strong phase arising from the radiative corrections of 0(a s ). It has been shown 
that using 5' = 30° together with the parameters given in Eq. (|28jl . one can explain the large 
negative value of sin(20i)0x s as well. Notice that the new coupling terms 0^22 an d ^222 are 
relevant only to the process b — > sss, so they do not affect other B — > PP and VP decays, 
such as B — > 7T7T, 7rif , pvr, pi^, etc, which are already well understood within the SM. 

In the case of the large Xa,h effects with Xmm = 7-5 where the BR of B + — > 77'i^ + is large, 
we can use the R-parity violating SUSY couplings to raise the BR of B + — > (which is 
small, 4.02 x 1(T 6 to begin with). It is possible to raise B(B + -> 0K+) to (8 - 9) x 10" 6 . 
However, in this case, sin (2(pi), Ks can not be large negative j^. 

The RPV terms can arise in the context of SO (10) models which explain the small 
neutrino mass and has an intermediate breaking scale where B — L symmetry gets broken 
by (16 + 16) Higgs. These additional Higgs form operators like 16#16 m 16 m 16 m /M p ; (16 m 
contains matter fields) and generate the RPV terms 

(b) R-parity conserving SUSY case 

As an example of the RPC SUSY case, we will consider the supergravity (SUGRA) model 
with the simplest possible non-universal soft terms which is the simplest extension of the 
minimal SUGRA (mSUGRA) model. In this model the lightest SUSY particle is stable and 
this particle can explain the dark matter content of the universe. The recent WMAP result 
provides [29]: 

ficoM/i 2 = 0.11261°;°°*, (29) 

and we implement 2a bound in our calculation. 

In the SUGRA model, the superpotential and soft SUSY breaking terms at the grand 
unified theory (GUT) scale are given by 



Y U QH 2 U + Y D QH 1 D + Y L LH 1 E + /iH^, 
17 



£ so ft = -^ m il^| 2_ \^m a KK + BiiH x H 2 

+(A U QH 2 U + A D QH l D + A L LH 1 E) + H.c] , (30) 

where E, U and D are respectively the lepton, up-quark and down-quark singlet superfields, 
L and Q are the SU(2)l doublet lepton and quark superfields, and H\p are the Higgs dou- 
blets, (pi and X a denote all the scalar fields and gaugino fields, respectively. In the mSUGRA 
model, a universal scalar mass m , a universal gaugino mass m-1/2, and the universal trilinear 
coupling A terms are introduced at the GUT scale: 

m\ = ml, m a = m 1/2 , A U ' D ' L = A Y U ' D ' L } (31) 

where Y U,D,L are the diagonalized 3x3 Yukawa matrices. In this model, there are four 
free parameters, m , my2, A , and tan/3 = (iT 2 }/(-ffi), in addition to the sign of fi. The 
parameters m^, // and A can be complex, and four phases appear: 9 a (from A ), 81 (from 
the gaugino mass mi), #3 (from the gaugino mass 7713), and 9^ (from the \x term). 

It has been shown in Refs. jsoL that the mSUGRA model can not explain the large 
negative value of sin(20i)^ s , because in this model the only source of flavor violation is in 
the CKM matrix, which can not provide a sufficient amount of flavor violation needed for 
the b — ► s transition in the processes B — > cf)K. The minimal extension of the mSUGRA has 



been studied to solve the large negative sin(2^i)^ s in the context of QCDF |30J], or both 
large negative sin(20 1 ) < ^^ s and large BR of B — > rj'K in the context of NF [31^ . 

The minimal extension of the mSUGRA model contains non-universal soft breaking A 
terms, in addition to the parameters in the mSUGRA model. In order to enhance contribu- 
tions to the b — > s transition, the simplest choice is to consider only non-zero (2,3) elements 
in A terms which enhance the left-right mixing of the second and third generation. The A 
terms with only non-zero (2,3) elements can be expressed as 

A U,D = A Y U > D + AA U ' D , (32) 



U,D 



■jU,D 



e 1 ^ with 



U,D 





U,D _ 



where AA U,D are 3x3 complex matrices and AA U { - D 

unless = (2, 3) or (3, 2). It is obvious that the mSUGRA model is recovered if AA 
0. 

For our analysis, we consider all the known experimental constraints on the parameter 
space of the model, as in Ref. (3(|. Those constraints come from the radiative B decay 

18 



TABLE IV: The branching ratios (in units of 1(T 6 ) for B + -> rj'K + (left) and B + -> (right) 



at tan/3 = 10 with non-zero Aj4.23 and AA^. The units for m 1( / 2 , |A)|) an d 



1) 



^ /i 23(32) 



are in GeV. 



\M 


800 


600 


400 





£i/1 23(32) 


m i/2 = 300 


79.6 


9.9 


81.0 


9.2 


79.6 


9.1 


79.0 


8.1 


66-74 


mi/2 = 400 


78.2 


9.9 


83.0 


9.6 


79.0 


9.2 


81.0 


8.5 


150 - 168 


m i/2 = 500 


84.8 


9.9 


83.7 


9.9 


81.0 


10.0 


77.0 


8.1 


244 - 256 


m i/2 = 600 


73.0 


7.6 


71.0 


7.5 


70.0 


7.5 


70.0 


7.1 


270 - 304 



process B -> X s7 (2.2 x 10" 4 < B(B -> X s7 ) < 4.5 x 10" 



), neutron and electron 



electric dipole moments (d n < 6.3 x 10 26 e cm, d e < 0.21 x 10 26 e cm js^l), reuc density 
measurements, K°-K° mixing (AM K = (3.490±0.006) x 10" 12 MeV Q), LEP bounds on 
masses of SUSY particles and the lightest Higgs (m^ > 114 GeV). From the experimental 



constraints, we find that 9± « 22°, 9 3 « 30°, and 0. 



•11°. For the phase 9 a, we set 



9 a — t 1 "- K h as been noticed (30J that the SUSY contribution mainly affects the Wilson 
coefficients Cs s (7 7 ) and Cs g (7j) and these coefficients do not change the weak annihilation 
effects arising from the SM calculation. 

In our calculation, we consider the case with non-zero AA 23 and non-zero AA^ 2 for 
tan/? = 10. All the other elements in AA U ' D are set to be zero. We compute the BRs for 
B -> rj'K and B -> <pK in the case of AAf 3 ~ AAf 2 and </>f 3 7^ </>f 2 with tan/3 = 10. 
Table IV shows the BRs for B + — > 7/lT + and _B + — > calculated for various values of 
the parameters m-1/2, \Aq\ and |AA^/ 32 J. For each m\j% and \A§\, the left column shows the 
BR for B + — > r]'K + and the right column shows the BR for B + — > <pK + . All the predicted 
BRs in the table are well consistent with the experimental data. The BR for B + — > r/if" 1 " 
is estimated as (3.1 ~ 4.4) x 10 -6 which also agrees with the data. The higher tan/3 values 
are also allowed, but the al.owed range of m m becomes smal.e. We satisfy the rehc den 81ty 
constraint using the stau-neutralino co-annihilation channel 35J. 

For the numerical calculation, we used the QCD parameters given in Eq. (|T9~)) and the 
additional strong phase 5' = 0. The value of mu2 varies from 300 GeV to 600 GeV, and the 
value of \Aq\ varies from to 800 GeV. Even though the value of m is not explicitly shown, 
it is chosen for different m 1( / 2 and A such that the relic density constraint is satisfied, e.g., 
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for m 1 / 2 =300 GeV, mo varies in the range (70 — 110) GeV. The value of mo increases as 
my 2 increases. The value of AA^^) increases as my 2 does. The phases 0^ and 4>® 2 are 
approximately —40° to —15° and 165° to 180°, respectively. So far we have assumed that 
AA^ 3 32 = 0. But if we use Av4 2 / 3 32 7^ and AA^ 32 = 0, the value of sin(20)^ s is mostly 
positive. 

In passing, we note that the set of the same parameters used in our calculation can also 
produce the large negative value of sin(20i)^ s |27(. As a final comment, we note that in the 
case of the large Xa,h effects with Xmin = 7-5, it is possible to raise the BR for B + — > <fiK + 
to (8 — 9) x 10~ 6 . However, in that case, the large negative value of sin(20i)^ s can not be 
obtained fej 



V. CONCLUSION 

We investigated the decay processes B + ^ — ► rj'K + ^ in the QCDF approach. In or- 
der to reliably estimate the weak annihilation parameter Xa = (l + Pa^ a ^) m ^f" an d the 
hard spectator scattering parameter Xh = (l + Ph^ h ^) m ^ff" ar i sm S from logarithmic di- 
vergences in the end-point region, we used the global analysis for twelve B — > PP and VP 
decay modes, such as B — > ttti, ttK, pir, pK, urn, ujK. From the global analysis, we found 
that both the large effect of Xa,h (less reliable) and the small effect of Xa,h (more reli- 
able) are allowed. For the former case, the parameters pa,h and 4>a,h are determined to be: 
pY = PY = PY = PY = 1, YY = 5T°, YY = 52°, YY = -23°, YY = 180°. For the 
latter case, the parameters pa,h and 4>a,h are: p p A p = 0, p v A p = 0.5, = 1, p^ p = 0.746, 

YY = -6°, YY = YY = 180°. 

In the case of the large Xh,a effects allowed by the "good" fit (with Xmin = 8-6 for the 
twelve decay modes), the BRs for B + ^ — > r]'K + ^ and B + ^ — > (f)K + ^ calculated within 
the SM saturate the large values of the experimental results measured by Belle, BaBar, and 
CLEO. Thus, there is no room for invoking new physics effects on the quark- level process 
b — > sss, which are implied by the large negative value of sin(20i)^ s recently reported by 
Belle. 

In contrast, in the case of the small Xh,a effects that is theoretically more reliable, the SM 
prediction for these BRs is smaller than the experimental data. Since both B + ^ — > ri'K + ^ 
and -> 4>K + ^ have the same (dominant) internal process b — ► sss, we took into 
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account possible new physics effects on the b — > sss transition, as in 0,|27| for explaining 
the recent Belle measurement of sin(20i)^ s . Specifically, we considered two new physics 
scenarios: R-parity violating SUSY and R-parity conserving SUSY. In the RPV SUSY 
case, the BRs for -> rj'K+W are predicted as 73.9(67.8) x 1(T 6 and the BRs for 

B+(°) -> 0K+(°) are 10.2(9.5) x 10~ 6 which are consistent with the data. The relevant new 
couplings are found to be: |A 322 | = 0.086, |A 332 | = 0.089, |A 323 | = 0.030, L = 0.66, R = 
—2.25. As an example of the RPC SUSY case, we adopted the simplest extension of the 
mSUGRA model, which contains only non-zero (2,3) elements in the soft breaking trilinear 
coupling A terms, in addition to the other parameters of the mSUGRA model. Considering 
all the known constraints on the relevant parameter space, we found that for tan/3 = 10, 
B(B + -> r)'K+) = (70.0 ~ 84.8) x 10~ 6 and B(B + -> <pK + ) = (7.1 ~ 10.0) x 10" 6 , which 
are in good agreement with the experimental data. 



ACKNOWLEDGEMENTS 

The work of B.D was supported by Natural Sciences and Engineering Research Council of 
Canada. The work of C.S.K. was supported in part by Grant No. R02-2003-000-10050-0 
from BRP of the KOSEF and in part by CHEP-SRC Program. The work of S.O. and G.Z. 
was supported by the Japan Society for the Promotion of Science (JSPS). 



[1] A. Gordon et al. [Belle Collaboration], Phys. Lett. B 542, 183 (2002); T. Tomura [Belle 
Collaboration], talk at the 38th Rencontres de Moriond: Electroweak Interactions and Unified 
Theories, Les Arcs, France, March 15-22, 2003, to appear in the proceedings; A. Kuzmin [Belle 
Collaboration], talk at the 38th Rencontres de Moriond: QCD and Hadronic Interactions, Les 
Arcs, France, March 22-29, 2003, to appear in the proceedings. 

[2] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 91, 161801 (2003); 
|arXiv:riep-ex/0309025| |arXiv:hep-ex/0311016| M. Pivk [BABAR Collaboration], talk at the 
38th Rencontres de Moriond: QCD and Hadronic Interactions, Les Arcs, France, March 22-29, 
2003, to appear in the proceedings. 

[3] D. Cronin-Hennessy et al. [CLEO Collaboration], Phys. Rev. Lett. 85, 515 (2000); 



21 



S. J. Richichi et al. [CLEO Collaboration], ibid. 85, 520 (2000); C. P. Jessop et al. [CLEO 
Collaboration], ibid. 85, 2881 (2000); R. A. Briere et al. [CLEO Collaboration], ibid. 86, 3718 
(2001). 

[4] D. Atwood and A. Soni, Phys. Lett. B 405, 150 (1997). 
[5] A. L. Kagan and A. A. Petrov, arXiv:hep-ph/9707354 



[6] M. R. Ahmady, E. Kou and A. Sugamoto, Phys. Rev. D 58, 014015 (1998). 

[7] For an application of the anomalous g — g — rj coupling to A& decays, please look at 
M. R. Ahmady, C. S. Kim, S. Oh and C. Yu, |arXiv:hep-ph/030503l"| 

[8] I. Halperin and A. Zhitnitsky, Phys. Rev. Lett. 80, 438 (1998). 

[9] A. Ali, J. Chay, C. Greub and P. Ko, Phys. Lett. B 424, 161 (1998). 

[10] N. G. Deshpande, B. Dutta and S. Oh, Phys. Rev. D 57, 5723 (1998). 

[11] D. S. Du, C. S. Kim and Y. D. Yang, Phys. Lett. B 426, 133 (1998). 

[12] M.-Z. Yang and Y.-D. Yang, Nucl. Phys. B 609, 469 (2001). 

[13] M. Beneke and M. Neubert, Nucl. Phys. B 651, 225 (2003). 

[14] E. Kou and A. I. Sanda, Phys. Lett. B 525, 240 (2002). 

[15] D. Choudhury, B. Dutta and A. Kundu, Phys. Lett. B 456, 185 (1999). 

[16] B. Dutta, C. S. Kim and S. Oh, Phys. Lett. B 535, 249 (2002). 

[17] B. Dutta, C. S. Kim and S. Oh, Phys. Rev. Lett. 90, 011801 (2003). 



[18] A. Kundu and T. Mitra, arXiv:hep-ph/0302123 



[19] Z.-j. Xiao, K.-T. Chao and C. S. Li, Phys. Rev. D 65, 114021 (2002). 

[20] A. Ali and C. Greub, Phys. Rev. D 57, 2996 (1998); Y.-H. Chen, H.-Y. Cheng, B. Tseng, and 

K.-C. Yang, Phys. Rev. D 60, 094014 (1999). 
[21] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nucl. Phys. B 606, 245 (2001). 
[22] D. Du, J. Sun, D. Yang and G. Zhu, Phys. Rev. D 67, 014023 (2003). 



[23] K. Abe et al. [Belle Collaboration], |arXiv:hep -ex/0308035 

[24] D. Du, H. Gong, J. Sun, D. Yang and G. Zhu, Phys. Rev. D 65, 074001 (2002); ibid. 65, 

094025 (2002) [Erratum-ibid. D 66, 079904 (2002)]. 
[25] T. Feldmann, P. Kroll and B. Stech, Phys. Rev. D 58, 114006 (1998). 
[26] M. Beneke and M. Neubert, Nucl. Phys. B 675, 333 (2003). 
[27] B. Dutta, C. S. Kim, S. Oh and G. Zhu, |arXiv:hep-ph/0312389l 



[28] R. N. Mohapatra, arXiv:hep-ph/9604414, 

22 



[29] C. L. Bennett et al, |arXiv:astro -ph/0302207 . 

[30] R. Arnowitt, B. Dutta and B. Hu, Phys. Rev. D 68, 075008 (2003). 

[31] S. Khalil and E. Kou, Phys. Rev. D 67, 055009 (2003); |a7^v:hep-ph/0307024l 

[32] M. S. Alam et al. [CLEO Collaboration], Phys. Rev. Lett. 74, 2885 (1995). 

[33] A. J. Buras, A. Czarnecki, M. Misiak and J. Urban, Nucl. Phys. B 631, 219 (2002). 

[34] K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002). 

[35] J. R. Ellis, T. Falk, G. Ganis, K. A. Olive and M. Srednicki, Phys. Lett. B 510, 236 (2001); 
R. Arnowitt, B. Dutta and Y. Santoso, Nucl. Phys. B 606, 59 (2001); M. E. Gomez and 
J. D. Vergados, Phys. Lett. B 512, 252 (2001); H. Baer, C. Balazs, A. Belyaev, J. K. Mizukoshi, 
X. Tata and Y. Wang, arXiv:hep-ph/0210441 ; A. B. Lahanas and D. V. Nanopoulos, 
arXiv:hep-ph/0303130"l 



23 



